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Effective field theory of BFKL pomerons interacting by QCD triple pomeron vertices is investigated. 
Classical equations of motion for the effective pomeron fields are presented being a minimal extension of the 

O-i! 

Balitsky-Kovchegov equation that incorporates both merging and splitting of the pomerons and that is self- 

| dual. The equations are solved for symmetric boundary conditions. The solutions provide the dominant 

contribution to the scattering amplitudes in the semi-classical approximation. We find that for rapidities 

of the scattering larger than a critical value Y c at least two classical solutions exist. Curiously, for each 

' of the two classical solutions with the lowest action the symmetry between the projectile and the target 

5— i ' 

is found to be spontaneously broken, being however preserved for the complete set of classical solutions. 
The solving configurations at rapidities Y > Y c consist of a Gribov field being strongly suppressed even at 
very large gluon momenta and the complementary Gribov field that converges at high Y to a solution of 
Balitsky-Kovchegov equation. Interpretation of the results is given and possible consequences are shortly 
discussed. 



'Email: sergb@mail.desy.de 
^E-mail: motyka@th.if.uj.edu.pl 



1 Introduction 



Understanding of scattering of hadrons and nuclei at very high energies in terms of Quantum Chromody- 
namics remains one of the most important challenges for the theory of strong interactions. The reasons for 
that are both of the phenomenological and of the purely theoretical nature. In practical terms, the basic 
physics of the present and future colliders, HERA, Tevatron, RHIC and the CERN Large Hadron Collider 
(LHC), is the physics of high energy scattering in QCD. On the other hand, the high energy regime of QCD 
reveals intriguing similarities to high energy regime of string theory pQ. To trace and explore postulated 
dualities between string theories and gauge theories and possible manifestations of those dualities in high 
energy scattering is a goal of primary importance. 

The most successful approach to high energy scattering in QCD is based on infinite resummations of large 
logarithms of collision energy ^/s in perturbative expansions of the scattering amplitudes. The corner stone 
of this formalism is the evolution equation derived by Balitsky, Fadin, Kuraev and Lipatov (BFKL) [2J |3J |I] . 
In the BFKL framework, the energy evolution of the exchange of an interacting gluon pair in a colour singlet 
state (the BFKL pomeron) was analyzed at the leading logarithmic (LL) approximation, and consequently, 
the energy dependence of the hard scattering amplitude was obtained that exhibits a power like growth 
with energy, A ~ s 1+A with the intercept A ~ 0.3. Such behaviour would eventually lead to violation of 
unitarity. Clearly, this indicates that when the energy is sufficiently large unitarity corrections to the BFKL 
evolution must be added. 

In recent years unitarity effects in high energy QCD have been vigorously investigated along two main 
lines. The Color Glass Condensate (CGC) approach Ej is formulated in the transverse position space 
and it is based on energy evolution of Wilson loops and phenomena of rescattering and recombination. In 
the large N c limit, the dynamics of the Color Glass Condensate may be analyzed in terms of a statistical 
model of color dipoles [7]. At a very general level, the dipole description of high energy scattering may 
be presented as a combination of multiple dipole splittings, the rescattering of dipoles off a target and a 
stochastic fluctuation term 

The QCD Reggeon Field Theory (QCD-RFT) approach is formulated in momentum space and bases the 
on standard diagrammatic calculus [§l 1101 ITT1 1121 1131 114j . The main building blocks here are QCD reggeon 
Green's functions and multi-reggeon vertices derived in the perturbative QCD. Scattering amplitudes may be 
represented in terms of Feynman diagrams of effective (non-local) reggeon fields. This effective field theory 
is constructed using the so-called Extended Generalized Leading Logarithmic Approximation (EGLLA) 
which resums the leading powers of a s log s for given topology of the reggeon diagram. The CGC and the 
QCD-RFT formulations should be, in fact, two different descriptions of the same theory, and they should 
be equivalent. 

In parallel to the theoretical efforts to determine the deep fundamental structure of the effective field 
theory for high energy scattering, more phenomenological studies of the unitarity effects have been carried 
out. Probably, a pair of the most fruitful (and entangled) concepts of the last decade were the Balitsky- 
Kovchegov (BK) evolution equation and the saturation model proposed by Golec-Biernat and 

Wiisthoff (GBW) [T3E1. The BK equation was derived in the CGC formulation and in the framework 
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a) b) c) 

Figure 1: Examples of diagrams of the effective field theory of QCD pomerons interacting with triple pomeron 
vertices: a) a fan diagram; b) a tree diagram defining the classical limit; c) a diagram with quantum loops. 

of QCD-RFT it can be viewed as a resummation of the BFKL pomeron fan diagrams of the type depicted 
in Fig. [JJi. An important feature of the BK equation is its simplicity - absence of the vertex for pomeron 
splitting removes all quantum loops of the complete theory. Therefore, a classical treatment of PFT is exact 
in the case of the BK equation. 

Numerous explicit solutions |191 1201 I21| and semi-analytical analyzes of the BK equation were pre- 
sented ^2 E21 [22] and applications were developed for the deep inelastic scattering (DIS) of nuclei and 
off the nucleon |241 125| I2fi| . Probably, the most remarkable features of the BK equation is generation of 
the saturation scale Q s (y), increasing exponentially with the available rapidity y: Q s {y) ~ exp(Ay), and 
the geometric scaling |27| I28[ I29j . Both the features of the BK equation are strongly favoured by the 
experimental data for o~{~f*(Q 2 )p) down to Q 2 = 0, the diffractive DIS data and heavy flavour electro- 
production ^JEHED]; two-photon cross sections [HJ; exclusive photo- and electroproduction of vector 
mesons and the deeply virtual Compton scattering jHSHHB]; etc. 

Unfortunately, the BK equation is not sufficient to address the problem of scattering of two similar 
objects, that is the symmetric situation, like proton-proton and nucleus-nucleus collisions at RHIC and the 
LHC. At present, efforts are made to understand the effective PFT beyond the BK limit [8ll3" H l35 ll36ll3T| l38]. 
The first key step is to determine properties of the theory after the vertex describing splitting of the pomeron 
to two pomerons is added to the BK framework. Importantly, the form of this vertex is imposed by the 
form of the vertex for two pomeron merging. The reason for that is that the two vertices differ only by 
the choice of direction of evolution in rapidity, which is arbitrary and should have no impact on physics. 
This simple fact was known in the RFT since a long time |391 140[ Hd] and recently it was re-discovered in 
the CGC framework, as the self-duality of the CGC Hamiltonian 02]- With the pomeron splitting vertex 
quantum pomeron loops become possible, as exemplified in Fig. ^i, and solving the theory becomes much 
more difficult. There exists, however, an interesting limit in which some explicit results can be found. 
Namely, collision of two large nuclei composed of A S> 1 nucleons each, can be analyzed. In the limit of 
very large A the quantum pomeron loops will provide only subleading contribution, suppressed by powers 
1/A. Then, the tree topologies of the pomeron diagrams give the dominant contribution to the 5-matrix, 
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corresponding to the classical limit of the effective Pomeron Field Theory. A pomeron diagram that has 
the tree topology is shown in Fig.^o. This concept was put forward and developed in a series of pioneering 
papers by Braun |34[ I35[ I36| , and an analogous path was taken by Balitsky in the framework of the Color 
Glass Condensate |38| . In the present study we provide a complementary analysis to that of Braun, solving 
the classical equations of motion of the Pomeron Field Theory, and reporting observations of some new, 
unexpected features of the PFT. Among them, perhaps, the most surprising is a break-down of the symmetry 
between the target and the projectile at the classical level. 

The paper is organized as follows. In the next section we describe the formalism of the effective QCD 
Pomeron Field Theory. In Sec. 3 we provide insight into the theory coming from a toy model of the Reggeon 
Field Theory in zero transverse dimensions. In Sec. 4 solutions of PFT are presented. We discuss the results 
in Sec. 5 and conclude in Sec. 6. 



2 Formalism 

2.1 The Balitsky-Kovchegov equation and the triple pomeron vertex 

Scattering of the small perturbative probe off a large nucleus was studied by Balitsky and Kovchegov Hoi 
116) . Balitsky derived the rapidity evolution equations describing the scattering amplitude in QCD, which 
involved an infinite tower of the correlators. This hierarchy is cut down to the lowest correlator in the large 
N c limit, leading to a much simpler and better tractable evolution equation. Such an equation was derived 
by Kovchegov in the framework of the dipole model. Thus, the Balitsky-Kovchegov equation provides the 
correct limit of the scattering amplitude in QCD for a large nucleus at the LLs accuracy, and for \/N c —* 0. 
In the infinite momentum frame, the BK equation resums the BFKL pomeron fan diagrams (see Fig. ^,), 
with the triple pomeron vertex equivalent to the Bartels vertex at the leading 1/N C accuracy. 

The BK equation was initially proposed as a non-linear evolution equation for the dipole scattering 
amplitude N(y;r, b), where the dipole spans the vector r and is located at the transverse position b. Thus, 
the BK equation reads 

dN(y;r,b) 



. hj a s (K®N)(y;r,b) 

where a s = N c a s /TT, and the linear term is determined by the BFKL kernel in the position space 

{It ® N)(y; r, b) = J ^ j N(y; r + r', b) - ~ 2 - £ + ^ N(y; r, b) | . (2) 

In the limit of the small scattering amplitude the quadratic term may be neglected and the BFKL equation 
in the dipole picture is obtained. 

The BK equation is a differentio-integral equation, for which the integral kernel depends on two two- 
dimensional vectors r and b. Numerical solution of the complete equation is possible [20] but cumbersome. 
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Besides that, the treatment of the large dipoles and large impact parameters in the BK equation is, strictly 
speaking, incorrect as far as QCD is concerned. Namely, the confinement of colour is not accounted for, 
which can be seen, for instance from the conformal invariance of the equation. Thus, the Froissart limit for 
the scattering matrix is broken due to a power like diffusion to the large impact parameters |43j . Thus, in 
this work we follow the approximation made in the earlier analyzes that the dominant contribution to the 
scattering amplitude comes from perturbative dipoles for which r <C R, where R is the nucleus size. In this 
case the evolution may be assumed to be local in the transverse plane and Eq. (^Q) becomes independent of 
b. Thus, the b-dependence can be suppressed in N(y; r, b) and it enters only through the initial condition 
for the evolution equation. 

For an azimuthally symmetric solution, N(y,r) = N(y,r), it is convenient to Fourier transform Eq. Q 
to the momentum space, 



<P(k 2 ,y) = I ^expHk.r)^%^ = f° - J (kr) N(y,r), (•">) 
■ 2tt r z Jn r 



where Jo is the Bessel function. In this case the following equation is obtained 

= a s (K'®<P)(y,k 2 ) - a s <P 2 (y,k 2 ), (4) 

and the action of the BFKL kernel (suitably shifted in the space of the Mellin moments in k 2 ) is given by 
(r>^A\( i*\ rda 2 j a 2 ^^ 2 ) - k 2 4>(y,k 2 ) k 2 <t>(y,k 2 ) \ 

where k 2 and a 2 are the virtualities of the exchanged gluons in the BFKL ladder. 

Equation Q may be further transformed to the form dependent on the unintegrated gluon density in 
the transverse space |25| . One has 



f(y,k 2 ) = -^k 4 V 2 cj ) (y,k 2 ), (6) 



and conversely, 



The unintegrated gluon density in the transverse space may be related in the small- a; limit to the collinear 
gluon distribution of the target A 

where y = log(l/x) and we restored the dependence of f(y, k 2 , b) on the tranverse position b, assuming that 
it is mild enough to ensure effective decoupling of the b dependence from the BK evolution, which should 
hold for a large target. 

After this transform is executed the BK equation reads 



d y f(y,k 2 ) = ^k 2 f ' hl 



7r J a 



2 



f(a 2 ) - f(k 2 ) + f(k 2 



a 



2 ~ k2 \ [4a 4 + fc 4 ]5 
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a) b) c) d) 



Figure 2: Elements of the effective action: a) the BFKL pomeron propagator; b) the merging vertex; c) the 
splitting vertex and d) the external sources of the fields. The arrow indicates the direction of evolution. 



2vra 2 



L % m ° 2) L £ b2) + ^ k2) L % ^ (i) ° 2) ] • (9) 

It is easy to verify that the nonlinear term describing joining of two pomerons, (/, /) — ► / is generated from 
the amplitude of the Bartels triple pomeron vertex (in the forward limit), 

(f\V 3P \f®f) = -2™* [ d 4a 2 fHv,a 2 ) f ^ f^) J \^ f{y^) 



- 2nal J f\y, a 2 )/( y , a 2 ) £ ^ log ^ j f(y, 6 2 ). (10) 

by functional differentiation with respect to the pomeron field P(y, k 2 ). The details of the complete effective 
action for interacting pomerons are given in the next section. 

2.2 Effective action and the self-duality 

Following Braun |34l 1351 05] . we shall construct an effective action to represent both pomeron merging and 
splitting. We wish to study the minimal extension of the BK equation, thus we neglect vertices at which 
more than three pomerons meet. It is clear that the vertex for pomeron splitting must be the same as the 
vertex for pomeron merging. The reason is that in order to distinguish merging of the pomerons from a 
splitting of a pomeron we have to specify the direction of the evolution in rapidity. This choice is, however, 
completely arbitrary and the form the action should not depend on it. More precisely, the inversion of the 
direction of evolution, y — ► — y results with the interchange of the outgoing (/t) and incoming (/) pomeron 
fields. In result, the effective action is invariant under the transform, 

f^f f , y^-y, (ii) 

and the form of the splitting vertex (p (g> P\V^ p \f) is given by (/^|Vap|/ ® f)^. 

Thus, in what follows we shall write the effective action for the interacting pomeron system based on 
the following principles: 
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1. Free propagation of the pomeron fields is described by the forward BFKL equation. 

2. Pomerons interact only through triple pomeron vertices. 

3. The vertex for merging of the pomerons is the Bartels triple pomeron vertex at large N c and in the 
forward limit, equivalent to the vertex that generates the Balitsky-Kovchegov equation. 

4. Splitting and merging of the pomerons are described by identical vertices. 

5. The fields do not depend on the transverse position. 

The elements of the action are graphically represented in Fig. |2j The above assumptions lead to a unique 
form of the action, 

A[f, /+; Y] = j^-j £ dy {C [f, /+] + C 3 [f, /t] + c\[f, /t] + C E [f, /t]} , (12) 
where the Lagrange function for the free propagation reads 



1 [da 2 



2 j a 

da 2 fdb 2 



f(y, a 2 )d y p(y, a 2 ) - f\y, a 2 )d y f(y, a 2 ) 



+ /^f J^fHy,a 2 )}Co(a 2 ,b 2 )f(y,b 2 ), 
and /Cq is the amputated forward BFKL kernel given by 



-rrJCo {a ,b )f{b = a 

4 7T 



b 2 



f(b 2 ) - f(a 2 ) f(a 2 ) 



\b 2 -a 2 



[46 4 + a 4 ]2 



The Lagrange function describing merging of two pomerons takes the form 

C 3 [f, /t] = -2,a 2 s I d 4 a 2 P(v, a 2 ) f %■ f(y, b 2 ) f f(y, c 2 ) 

J CL J a 2 Ja 2 C 



27TQ„ 



da 



db 2 



Splitting of a pomeron contributes with 



4[f,p] = -2W S J ^a 2 f(y,a 2 ) £ ^ f\y,b 2 ) £ * - f\y,c 2 

-2na 2 J^f(y,a 2 )fHy,a 2 ) £ & log (^) f\y,b 2 ). 
The coupling of the pomerons to the external sources is represented by 

c E [f,f*] = J [fky,a 2 )f(y,a 2 ) + fHy,a 2 )f E (y,a 2 ) , 

where the external sources will be assumed to be localized in rapidity, 

f E (y,a 2 ) = f A (a 2 )S(y), f E (y,a 2 ) = fg(a 2 )5(y — Y). 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 
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Clearly, /a represents the amplitude of emission of the pomeron described by the field f(y,k 2 ) from the 
source at y = and y B is the coupling of f^(y, k 2 ) to an external source at y = Y. 

Let us stress that the treatment of the transverse position is missing and thus the evaluation of the 
quantum loops of the complete theory cannot be performed with this action. Nevertheless, the action 
treated in the semi-classical framework may be used to approximately resum the BFKL pomeron tree 
diagrams (see Fig. ^p) in scattering of two large objects, for instance of two nuclei. For a scattering in which 
the projectile and the target have sizes much larger than the typical momenta in the QCD pomeron, the 
momentum transfer of the pomeron line originating from the external particles is bounded to be small by the 
form-factors of the sources and may be therefore neglected. In the diagrams without closed pomeron loops 
the constraint imposed on momentum transfer of the external lines propagates and extends to all pomeron 
lines. 

In fact, a complete action that properly represents the degrees of freedom corresponding to the momen- 
tum transfer (or equivalently to the transverse positions of the pomerons), that are missing in the action 
(|12|) , was proposed by Braun [HE] • It is also straightforward to write down an analogous action in the present 
formulation. That complete action, however, leads to the same dynamics of a scattering of two large objects 
in the semi-classical limit, the problem that we address in this work. Therefore, we restrict ourselves to the 
simplified effective action given by (|12j) . 

Let us return to the symmetry of the action defined by Eq. (|11|) . This symmetry causes the action to 
be self-dual. Indeed, after integration by parts of the "time derivative" part of the action 

Y dyl\f(y,a 2 )d y f\ y ,a 2 ) - p(y,a 2 )d y f(y,a 2 )} - F dy \-f\y, a 2 )d y f(y, a 2 )} + (. . .) (19) 
2 I J Jo I J 

where (...) denote the boundary terms, one gets that 

This means, that the field P{y,k 2 ) is the canonical conjugate of f(y,k 2 ), up to the factor of l/k A which 
can be easily absorbed into the field definitions and trivial complex phase factors. After invoking the 
symmetry (|11[1 we conclude that the bulk part of the action (|12j) may be rewritten in the self-dual form. 
The symmetry of the action 1)12|) may be completed by assuming the symmetric external sources, that enter 
Ce- Then, one expects the solution of the field equations {/,/'} to be also symmetric 

f(y,k 2 ) = f\Y-y,k 2 ). (21) 

In what follows, we shall refer to this as to the projectile-target symmetry. 

2.3 Equations of motion 

Let us list the functional derivatives of the action of the effective Pomeron Field Theory (|12j) with respect 
to P{y,k 2 ): 

tr„\f ftl II 

:f(y,k 2 ); (22) 



o 



S(d y p(y,k 2 )) 2fc 4 ' 
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S£o[f,f] __ 1 
5p(y,k*) " k 4 



\d y f{y,k 2 )-j ^IC (k 2 ,b 2 )f(y,b 2 



(23) 



^ 3 [/,/t] 
6p(y,ki) 



; (24) 



^4[/./ f ] 

<5/t(y,fc 2 ) 



-2ira s p 



2f k2d 4^f(y^)f ^-f\y,b 2 ) 



&C g [/,/t] = 1 

sp(y,k 2 ) 



f E (y,k 2 



(25) 



(26) 



Analogously one computes the functional derivatives with respect to f(y,k 2 ). The results of that proce- 
dure may be obtained by an interchange of / f * in equations (|22H26j) and changing the sign of f(y, k 2 ) 
in (|22|) and of the d y f(y, k 2 ) in l|2H|) , Thus, one obtains the following equations of motion, 



d y f(y,k 2 ) = ^k 2 ft '" 2 



7T 



f(a 2 )-f(k 2 ) , f(k 2 ) 



+ 



-27raf 



-27ra c 



l« 2 -^ 2 l [4a 4 + fc 4 ]5 
da 2 



A* ij^ny^j/lfiy^ + f^k 2 )^ logf£l/ (y ,a 2 



2 jf £ a 2 /(, a 2 ) I £ IHV, + m ^) jf, £ log (| j ^ , 
-2vra 2 J Q k2 ^f(y,a 2 )f\ y ,a 2 ) log ^ , 



and 



5 fc 2 



7T 



da 2 
^2" 



/t(a 2 )-/t(fc2) , /t(fc 2 ) 



+ 



a 2 -/c 2 | [4a 4 + A; 4 ]5_ 



-2vtq' 



-27ra„ 



t2 L £ /,(! " ° 2) /, * ft2) + /tfe - * 2) /„ £ bs (f) - 2 > 



2 / ^L„ 2 ft^. .2 



<i& 2 



/(y,6 2 )+/t(y,A; 2 ) 



a 4 Jo? b 4 J k 2 a 4 



<ia 2 



log p /(y,a 2 ) 



2vra 2 J Q k2 ^f\y,a 2 )f(y,a 2 ) log ^ 



with the two-point boundary conditions, 

/(y = 0,fc 2 ) = / A (fc 2 ), 



/t( y = y,fe 2 ) = 4(fc 2 ). 



(27) 



(28) 



(29) 
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The equations are equivalent to the equations derived by Braun [31] , although they are formulated using 
other variables. Therefore we shall refer to equations (|27|). (|28(1 as to the Braun equations. Apparently, 
the present formulation is more complicated and less convenient than the original one. Still, it might be 
advantageous to use the present form. The reason is that the interpretation of the degrees of freedom that we 
use is straightforward in terms of perturbative QCD in the momentum space; the basic physical objects: the 
unintegrated gluon and the triple pomeron vertex in the momentum space are represented in a transparent 
way. Using the present form it should be also relatively simple to account for non-leading corrections to the 
BFKL kernel, as it was done for the BK equation |25j . 

2.4 The S-matrix 

Solutions to classical equations of motions for the pomeron fields may be used to determine the S-matrix for 
the high energy scattering in the semi-classical approximation. In order to do that, however, the dependence 
of the problem on the transverse position has to be taken into account. 

First, let us consider the general case, in which the action has the complete dependence on the transverse 
position. Thus, the pomeron fields depend on the position b\ and b<i with respect to the center of the 
projectile and the target, correspondingly: f(y,k 2 ) — > f(y,k 2 ,bi) and f\y,k 2 ) — > f^(y,k 2 ,b2). Suppose 
that we know the solutions to the Braun equations with the full dependence of the transverse position 
for a given impact parameter b of the collision. Then, the complete action may be evaluated for such a 
solution by performing integrations over transverse positions of all the fields, with the weights provided by 
the b-dependent Lagrangian density. 

In this paper we do not attempt to resolve the complex dynamics of the fields in the transverse plane. 
Therefore we should apply an approximate treatment, in the spirit of the original work of Kovchegov and 
following the initial Braun proposal. Those authors assumed that the sources of pomeron fields were large 
nuclei. Those objects are much larger than typical pomeron sizes, defined as inverse of the typical gluon 
virtualities in the pomeron. Therefore, for the bulk of interactions, an approximate translational invariance 
in the transverse space holds. This is not true only at the nucleus boundary, which gives, however, only a 
subleading contribution to the scattering amplitude. The simplest way to approximately account for this 
situation is to assume that the action is local in the transverse position. Then, it is enough to solve the Braun 
equations with input conditions dependent on the transverse position, /^(A; 2 , b\) and fg(k 2 , b\ — b) at given 
impact parameter vector b. In writing so, we assume that the initial condition Ja (k 2 , b±) is centered at b\ = 
and the distribution fg(k 2 ,b\ — b) develops around the point whose position is given by b. For instance, 
for a collision of two cylindrical nuclei with the same radius R, one has: /^(/c 2 , b\) = fA(k 2 )<d(R — b%) and 



Thus, assuming locality of the evolution in the transverse space the complete action takes the form 



P B (k 2 ,b 2 ) = fl(k 2 )Q(R-b 2 ). 



A[f,P;Y,b] = J d 2 b 1 A[f(y,k 2 ,b 1 )J~\y,k 2 ,b-b 1 y,Y] 



(30) 



where the equations of motion may be employed to obtain 



A[f, f; Y\ = \ J o _ dy {c E [f, f] - C 3 [f, f] - 4[f, f]} 



(31) 
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leading, in the semi-classical approximation, to the S"-matrix 

S(Y,b) = e W {-A[f,P;Y,b}} 



(32) 



Some more details and subtleties of this approximation will be discussed in the next section. 

3 Toy model — Reggeon Field Theory in zero transverse dimensions 
3.1 Formulation 

As a constructive example of a possible qualitative behavior of the interacting pomeron system we consider a 
zero dimensional toy model of interacting pomerons, the, so called, Reggeon Field Theory in zero transverse 
dimensions (RFT-0). The model of RFT-0 was formulated and studied in depth long time ago, see e.g. 
|39l !4Uj and recently it has enjoyed a revived interest |441 145| . In fact, it turns out that RFT-0 in the weak 
coupling regime exhibits some generic features which seem to be present also for interacting QCD pomerons. 
Therefore, this much simpler model may be used to provide some insight into the complex dynamics of QCD 
Pomeron Field Theory. 

This model is determined by the action 

A RF T-o[q{y),p{y)]Y} = [ dy £rft— o- (33) 

Jo 

with the Lagrangian: 

£rft-o = ^qdyp - ^pd y q + fiqp - Xq{q + p)p + p(y)q (y) + p (y) q(y) , (34) 

where \x is the intercept of the pomeron, A is the triple pomeron coupling and {q,p} are (up to complex 
phase factors) Gribov fields depending only on rapidity and responsible for the creation and annihilations 
of pomerons. They correspond to / and f* of the BFKL Pomeron Field Theory. The functions qo(y) 
and po(y) are the external sources of the q and p fields respectively. In analogy to the assumptions of the 
previous section we shall consider a scattering process at rapidity Y with the source terms assumed to take 
the form: 

Qo{y) = 01%), Po(y) = 92S(y-Y). (35) 
Note, that the action is invariant under the duality transformation 

p «-> q and y — > Y — y (36) 

for symmetric boundary conditions g% = 52 (obviously, the bulk action is invariant for any external cou- 
plings). 

Dynamics of the system defined by the Lagrangian (|34[) was intensively investigated both in the complete 
quantum framework [101 and in the semi-classical approximation jSH]- Here, we focus on the latter treatment 
in order to match the approximations which we use in the description of the BFKL Pomeron Field Theory. 
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Thus, the quantum evolution of the system may be represented by the ^-matrix expressed using the path 
integral 

S(Y;g 1 ,g 2 )= j [Dq Dp] exp {-A RFT - [q(y),p(y);Y}} , (37) 

where the probe trajectories obey q(0) = and p(Y) = 0, as the initial conditions are absorbed in the 
action. In the semi-classical limit, the dominant contribution to the path integral comes form the classical 
trajectories {q a ,p a } for which the action is stationary, 

S(Y; gi ,g 2 ) ~ ^2A a exp{-ARFT-o[q a ,P a ;Y]}, (38) 

a 

where A Q represent the quantum weights of subsequent classical trajectories, which at the leading approx- 
imation come from resummation of Gaussian quantum fluctuations around the classical trajectory. Note, 
that the system evolves in rapidity which is formally equivalent to an evolution in the Euclidean time, thus 
the 5-matrix is dominated by classical trajectories with the minimal value of the action. In is important 
to stress that the action of RFT-0 exhibits the feature of self-duality (or projectile-target symmetry) as its 
PFT counterpart. 

3.2 Solutions: spontaneous breaking of projectile-target symmetry 

The extremal value of the action is reached for classical trajectories {q,p} which obey the equations of 
motion, 

d y q = fiq — Xq 2 — 2Xqp (39) 
— d y p = up — Xp 2 — 2\qp (40) 

with the two-side boundary condition 

g(0) = gi , p(Y) = g 2 . (41) 

For <7i < fi/X and g 2 < p/X the classical trajectories are confined to a triangle in the phase space spanned 
by points with (p,q) coordinates: (0,0), (0,/i/A) and (fi/X,0). These boundary conditions permit for 
existence of multiple solutions provided that rapidity Y is large enough. Thus, for Y smaller than a 
critical value Y c (depending on g±, g 2 , X and /i) there exists a unique solution to the classical problem, 
{Qi(y'i <7i> 92),Pi(y', gi, 92)}- In the case of g\ = g 2 = g the solution preserves the symmetry between the 
target and the projectile, 

Qi(v,9,9) =Pi(Y-y;g,g). (42) 

This simple picture changes at Y = Y c . In this point two more solutions {q 2 (y; g, g),p 2 (y, g, g)} and 
{^2(2/; 9-> 9) > 2*2(2/; 9-> 9)} become possible which do not inherit the symmetry between the target and the 
projectile embedded in the action and the boundary conditions, 

Q2(y,9,g) ¥"P2{Y -y;g,g), and q' 2 (y,g,g) ^ p' 2 {Y - y;g,g). (43) 
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Figure 3: Classical solutions of the RFT-O: a) the {q,p} trajectories for Y > Y c ; b) value of the action 
ARFT-o[q(y}g,g),P(y)9>g)]Y] for the symmetric solution (dotted line) and the asymmetric solution (dashed 
line) as a function of scaled rapidity fjY. 

An example of the solutions is given in Fig. |3K plotted in the phase space {p, q} . The parameters of the 
model were chosen to be \xj\ = 5, g\ = <?2 = 0.7, and [iY = 8. For this rapidity, one finds the symmetric 
trajectory 1 and two asymmetric trajectories: 2 and 2'. At yet larger values of rapidity Y more solutions 
are possible, corresponding to cycles in the phase space and giving larger values of the action, so we neglect 
those cycles in the present analysis. 

The value of the action corresponding to trajectories 1 and 2 is plotted in Fig. Ob as a function of the total 
rescaled rapidity jiY . Note, that trajectory 2 is only possible for Y > Y c , and the critical rapidity Y c ~ 4 for 
our choice of parameters. Clearly, the value of the action is smaller for the asymmetric trajectories, therefore 
the asymmetric trajectories are expected to dominate the Euclidean path integral defining the scattering 
amplitude at large rapidities. At Y = Y c , however, the action of the asymmetric trajectory joins smoothly 
the action of the symmetric trajectory. Thus, one concludes that at the transition region of Y — Y c the 
contribution of trajectory 1 to the scattering amplitude should be also included. 

Note, that the emergence of the dominant asymmetric solutions may be interpreted in terms of sponta- 
neous breaking of a discrete symmetry of the action. The symmetry between the projectile and the target 
(leading to the self-duality of the action) is built in the action (|33|) and in the boundary conditions. Clearly, 
this is a discrete symmetry. The dominant solutions of the equations of motion, however, are not symmetric. 
Thus, the symmetry is spontaneously broken. This is possible, as the boundary conditions are defined at 
two points of rapidity and the classical solutions need not be unique. As usual, however, the symmetry still 
holds for the full set of solutions, 

(h{v,9,9) =P2{Y -y,g,g) and p 2 (y;g,g) = q 2 (Y - y;g,g). (44) 
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This means that the under the duality transformation (|H6|) each solution is transformed into itself (solution 1) 
or into another solution (solutions 2 and 2'), and the full set of solutions {p, q} is invariant under the duality 
transformation of p <-> q and y — > Y — y. 




Figure 4: Scattering amplitude iT = 1 — S from RFT-0 for g\ = g 2 = 0.7 at fi/X = 5 as a function of 
rescaled rapidity pY in various approximations: dotted line represents amplitude obtained from the symmetric 
solution 1, dashed line correspond to the pair of asymmetric solutions {2,2'}, the dash-dotted line accounts 
for sum of contributions from {1,2, 2'} and the solid line represents the full quantum solution of the problem, 
(the plot is taken from [461). 

The observed phenomenon of spontaneous symmetry breaking occurs at the classical level. At the 
quantum level, however, the projectile-target symmetry should hold 1 . One sees it, for instance, from the 
form of the 5-matrix in the semi-classical approximation using the three solutions {1, 2, 2'}. The calculation 
of the quantum weights for Y > Y c was performed in |3J leading to, 

S{Y;g 1 ,g 2 ) ~ - exp {-Arft-o[<7i,Pi; Y]} + exp {-Arft-o[Q2,P2; Y}} + exp {-Arft-oW 2 ,P2> Y }} > ( 45 ) 

where the minus sign of the first term comes form the complex phase factors picked up by the trajectory 1 

at the turning points. One sees that the symmetry between the projectile and the projectile is restored for 

the S-matrix already at the semi-classical level, by summation over the complete set of (asymmetric and 

symmetric) classical trajectories. It is interesting to ask, however, whether there could be some signs of the 

symmetry breaking found at the classical level, that would be seen for more exclusive observables, like for 

the rapidity distribution of produced particles. In principle, such a classical measurement should destroy 

the quantum coherence and select just one of the classical solutions. If this were true, it should lead to an 

asymmetric particle production between the identical projectile and target in individual events. 

x Due to quantum coherence in finite quantum systems the spontaneous symmetry breaking does not happen. It is only 
possible in the thermodynamical limit when the coherence between the asymmetric configurations is broken. 
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Figure 5: Comparison of the classical asymmetric solution of the full theory to the solution of the equation 
resumming the fan diagrams: a) the trajectory: q^iu) (upper solid line) and p2(y) (lower solid line) and 
the trajectories of the fan equation q~f{y) and pf(y) (upper and lower dotted lines respectively) for the total 
rapidity = 8 as a function of fiy; b) the action for the classical solutions: the asymmetric one (the 
dashed line) the symmetric one (the dotted line) versus the action resumming pomeron fan diagrams as a 
function of fiY . 



It is instructive to compare various approximations to the scattering amplitude iT = 1 — S. Thus, we 
plot in Fig. Qp, the scattering amplitude evaluated in the semi-classical approximation assuming that it is 
dominated by the symmetric solution 1 (which should be valid for Y < Y c ), by the pair of asymmetric 
solutions {2,2'}, (which should hold for Y 3> Y c ), and including contributions to the 5-matrix of all three 
solutions for Y > Y c . We expect the last and the most complete choice be the most accurate. The comparison 
to the exact quantum evaluation of the scattering amplitude 2 (the continuous curve) to the various semi- 
classical evaluations reveals, however, that the contribution of {2, 2'} approximates the exact answer best. 
It may be a coincidence or a hint on the subtle point of how to treat contributions to the S-matrix from 
subleading trajectories. We leave the issue for further studies. 



3.3 Fan dominance 



Anticipating the results of the next section, let us briefly mention an interesting feature of asymmetric 

solutions for Y being significantly larger than the critical rapidity Y c and for sources 51,32 <C /-t/A. Then, 

one of the fields, say q, grows with the increasing rapidity. The other field, p, assumed to take the value of 

92 *C jti/A at rapidity Y decreases further for decreasing y. Therefore the term involving p 2 in Lagrangian 

2 The curve was obtained in the course of our ongoing study of the RFT-0 |46|. where we solve the RFT-0 both at the classical 
and the quantum level. We leave the description of the details to that paper. 
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(|.'S4(1 is of little importance and it may be neglected. Thus, one ends up with the evolution equations of the 
system with the absent vertex for the pomeron splitting. This evolution equation resums the fan diagrams 
of the merging fields q, and it will be referred to as the fan equation. 

We exemplify in Fig. EH the similarity of the asymmetric solution of the full theory {<fe(y) , P2(y)} to the 
solution of the equation resumming the fan diagrams {qf(y),pf(y)}- The figure was obtained with fi/X = 5, 
9i = 92 = 0.7, and we set fiY = 8 ~ 2fiY c . The smaller fields P2(y) and p/(y) turned out to overlap 
with high accuracy. We find some noticeable difference between ^(y) and qj{y) only at higher values of 
fiy. We checked that the asymmetric solution is closer to the "fan equation" solution if Y is larger and gi 
and 52 are smaller. The results of evaluation of the action corresponding to {q~2(y) , P2(y)} , {Qf(y),Pf(y)} 
and to the symmetric classical solution {qi(y),pi(y)} are shown in Fig. 03. Again, the solution to the fan 
equation yields the action that reasonably well approximates the action Arft-o(<?2(?/; 9, g),P2(y, 9, g)', Y) of 
the dominant trajectory. 

We find this convergence of the classical system to the fan-dominated regime to be a rather surprising 
effect especially in view of the fact that we started from completely symmetric boundary conditions. A 
similar observation, however, was made in sligtly different realisation of RFT-0, with a non-zero four pomeron 
coupling 45 j. It will be very interesting to check whether a similar phenomenon occurs for the theory of 
interacting QCD pomerons. 



4 Solutions of the Braun equations 
4.1 Parameters and the solving procedure 

Classical equations of motion (|27|) and (|28|) for the effective pomeron fields f(y, k 2 ) and f^(y, k 2 ) were solved 
numerically for various values of the total rapidity Y of the scattering. We chose a fixed coupling constant 
a s = 0.2. The boundary conditions were assumed to be symmetric 

j. N k 4 

The form of the input condition was inspired by the properties of a saturated gluon distribution in the 
nucleon. Thus we set R 2 = 8 GeV -2 in order to match the preferred value of the nucleon size. The scale 
Qq = 0.5 GeV 2 corresponds to the saturation scale for 10~ 3 < x < 10~ 2 . The overall normalisation factor 
N = 2 so that the collinear gluon distribution obtained from the input is similar to the actual collinear 
gluon distribution xg(x, Q 2 ) in the proton for < x < and for moderate Q 2 . Let us stress that we 
made that choice only to pin down the physically relevant ranges of parameters. 

The numerical solution was based on a Chebyshev interpolation method in the variable log(/c 2 ) used to 
discretize the differentio- integral equations ()27|) and (|28|) . In order to solve the two-point boundary problem, 
we applied the iterative procedure defined by Braun [HSJ. Thus, in each iteration the evolution in rapidity of 
one of the field f(y, k 2 ) or f\y, k 2 ) was performed while the other field was set to its value obtained in the 
previous iteration. In the odd iterations, the field f(y, k 2 ) was evolved from its initial value /a(^ 2 ) from y = 
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to y = Y and the field f\y, k 2 ) was kept fixed (in the first iteration p(y, k 2 ) was set identically to zero). In 
the even iterations f^(y, k 2 ) was evolved from y = Y down to y = with the initial value fg{k 2 ) at y = Y. 
The iterations were continued until a fixed point of f(y, k 2 ) and P(y, k 2 ) was reached. Clearly, at the fixed 
point of the iterative procedure the fields {f(y, k 2 ), P{y, k 2 )} solve the system of Braun equations (|2*T|) and 
(|28|) with boundary conditions (|29|) . The method was found to be stable and robust and no problems with 
convergence occurred of the kind reported in [HE]- A disadvantage of the iterative method is that it may be 
only used to find the solutions which represent the attractive fixed points of the procedure. Unfortunately, 
we expect from the analysis of solutions of the RFT in zero transverse dimensions, that the solution to 
Braun equations is not unique at larger rapidities and there should exist multiple solutions. Therefore it is 
probable that the iterative method finds only some of them. On the other hand, one hopes that the most 
relevant solutions that minimize the action may become an attractive fixed points of a reasonable iterative 
procedure. Some more arguments in favour of this scenario will be given based on the properties of the 
found solutions. 

4.2 Properties of solutions and spontaneous symmetry breaking 

In what follows we shall describe the properties of the solutions to the Braun equations for rapidities of the 
scattering Y = 6,8, 10, 12 and Y = 16. The pomeron field f^(y,k 2 ) will be often presented as a function 
of the transformed rapidity y' = Y — y, so that the initial condition for f> is imposed at y' = 0. In the 
figures we shall plot f(y,k 2 )/k 2 and f^(y',k 2 )/k 2 unless we explicitly specify differently. Such a choice of 
variables is preferred by their solitonic behaviour in the case of the BK equation. For future reference, let 
us introduce the notation f BK (y,k 2 ) to represent the solution to the BK equation with the input given by 
(|46j). In the figures the label "Input" is used for fA{k 2 )/k 2 , where fA(k 2 ) defined by (HHJ). 

In Fig. |HK and Fig. |Eb we show the solutions to the Braun equations for Y = 6 and Y = 8 respectively. 
Solid lines denote f(y,k 2 )/k 2 and f\y' ,k 2 ) /k 2 is shown with points. The curves are plotted for y varying 
from zero to Y in the steps of one. Clearly, in both cases the solutions are symmetric, f(y, k 2 ) = f'(Y—y, k 2 ). 
Anticipating Fig. 112b . we point out that the the solutions exhibit a similar behaviour to solutions of the 
BK equation at large gluon momenta k 2 . At small momenta, below the saturation scale of the BK equation, 
f(y,k 2 ) and P(y', k 2 ) are much flatter than f BK (y, k 2 ). 

We find that the symmetry between / and /t breaks down at some critical rapidity Y c ~ 9. For Y > Y c 
only the asymmetric solutions are found, for which f(y,k 2 ) P(Y — y,k 2 ), compare Fig. EH (Fig. Eb) 
and Fig. El (Fig. EJ) for Y = 10 (Y = 16). Thus, the symmetry between the projectile and the target is 
spontaneously broken for individual classical solutions, in close analogy with the phenomenon appearing in 
RFT-0, described in detail in Sec. 13.21 The asymmetry between f(y, k 2 ) and p{Y—y, k 2 ) vanishes at Y = Y c , 
so the asymmetric solutions connect smoothly to the symmetric one at Y = Y c and the asymmetry builds 
up gradually with increasing Y. Certainly, the numeric value of the critical rapidity Y c is not universal, it 
depends on the boundary conditions and on the value of a s . It is important to note, that for each asymmetric 
solution {/, /t} there exists a complementary solution {/', /^'}, such that f'(y,k 2 ) = f^(Y — y,k 2 ) and 
/t (Y — y,k 2 ) = f(y,k 2 ), reflecting the symmetry between the projectile and the target encoded in the 
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action and the symmetric boundary conditions. Knowing that, in the further analysis of the solutions we 
choose arbitrarily that f(y,k 2 ) is the larger field and f^(y',k 2 ) is the smaller one. 

For Y > Y c , the general features of the larger field / are the following. At Y ~ Y c the solution is similar 
to the symmetric solutions found for Y < Y c . With increasing Y a pattern appears of a traveling wave, that 
is formation of a peak of f(y, k 2 )/k 2 traveling towards larger values of log(/c 2 ) with increasing rapidity with 
only small changes of the shape, see Fig. and Fig. 03. Recall, that it is behaviour characteristic for the 
BK equation |29j . The similarity of the solution to the BK solution will be investigated in more detail in 
Sec. PI 

The smaller field f*{y, k 2 ) evolves differently, see Fig.|Hland Fig. El At Y = Y c it matches f(Y — y,k 2 ) and 
for the increasing Y it experiences a significant overall suppression, stronger at larger Y. For instance, for 
y ~ Y/2 the maximal value of fi(y, k 2 )/k 2 is about an order of magnitude smaller than the maximal value of 
f(y,k 2 )/k 2 at Y = 10 (compare Fig. an d Fig. EJ) and about three orders of magnitude smaller at Y = 16 
(see Fig. 03 and Fig. EJ. Thus, we conjecture that in the limiting case of very large total rapidity Y, f* may 
is arbitrarily small except of the rapidities y ~ Y where the source term for f* is still important. In this 
context, it is instructive to study the y-dependence of P(y, k 2 )/k 2 at fixed k and compare it to f(y, k 2 )/k 2 . 
This comparison may be performed using Fig. I1UI It turns out, that there appear two distinct regimes of 
evolution of f^(y,k 2 )/k 2 with rapidity (at fixed momentum). Thus, if the field f(y,k 2 )/k 2 is strong, the 
smaller field f\y' ,k 2 )/k 2 is exponentially suppressed with increasing y', /'(y', k 2 )/k 2 ~ exp(— j3\y') with 
/?i ~ 1, crudely. This is the region where the absorption of f* by / drives the evolution of /'. Then, 
when y is sufficiently small and the field f(y,k 2 )/k 2 is weaker, the absorption becomes less relevant and 
P (y 1 ,k 2 ) jk 2 grows exponentially with increasing y', f^(y',k 2 )/k 2 ~ exp(/32y') with the exponent 02 — 0.4 
(with our choice of parameters) a value somewhat smaller than the BFKL intercept ujq = Aa s log(2) ~ 0.53. 
Note, that the characteristic rapidity y, at which the transition occurs from the strong absorption regime 
to the BFKL driven growth regime, depends on k. This is natural, as the field f(y,k 2 )/k 2 becomes strong 
at larger values of y for larger k. 

The shape of P{y' ,k 2 )/k 2 in k 2 exhibits some interesting features too, see Fig. Eland Fig- El At small 
values of k, P(y', k 2 )/k 2 tends to a flat function. This should be compared with the case of the BK where 
f BK (y,k 2 )/k 2 ~ k 2 at small k. On the other hand, at large k 2 the decrease of /'(y', k 2 )/k 2 with increasing 
k 2 is slower than the decrease of f BK (y,k 2 )/k 2 . Thus, the overall picture is that f^(y',k 2 )/k 2 is much 
flatter than / BK (y, k 2 )/k 2 . As rather surprising comes an observation that effects of non-linear interactions 
in f^(y',k 2 )/k 2 extend to very large values of momenta, causing a strong suppression of f^(y',k 2 )/k 2 for 
all momenta up to k = 10 3 GeV, the value larger than the saturation scale generated by the large field 
f(y,k 2 ), see for example Fig. El In fact, we checked that the suppression of f^(y',k 2 )/k 2 in comparison to 
f BK (y, k 2 )/k 2 is strong even at k = 10 5 GeV (not shown). 

The explanation of this phenomenon is the following. At large values of k 2 the input function fg(k 2 ) 
was assumed to tend to a constant, in other words the anomalous dimension vanished for the input. For the 
BFKL or the BK system the rapidity evolution generates an anomalous dimension of 70 — 0.3 — 0.5, strongly 
enhancing f(y, k 2 ) for large k 2 and y. For f*(y', k 2 ), however, the evolution and BFKL diffusion are almost 
completely blocked by large absorptive corrections coming from the interaction of p(y' \k 2 ) with the large 
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field f(y, k 2 ). Recall, that the input for f\y', k 2 ) resides in y = Y, where the value of the field f(y, k 2 ) is 
the largest and so is the related saturation scale. Therefore, before any BFKL diffusion or enhancement of 
P(y',k 2 ) becomes possible (that is at sufficiently small y') strong suppression of f^(y',k 2 ) occurs and the 
population of the large momenta region is initiated from a drastically reduced f'(y', k 2 ). 

In order to provide a more synthetic picture of the behaviour of f(y, k 2 ) and f\y, k 2 ) we illustrate the 
case of Y = 16 with three dimensional plots of the solutions shown in Fig. 111! Note, that we plot in this 
figure f^(y, k 2 ) instead of f\y' , k 2 ). Thus, the input for / appears at y = in Fig. II lb and the input of p 
is plotted for y = 16 in Fig. ITTb . 

4.3 BK fan dominance 

We have already related briefly the larger component f(y,k 2 ) of the solution to the Braun equation to 
the solution of the Balitsky-Kovchegov equation / BK (y, k 2 ). A more detailed comparison is performed for 
f(y,k 2 ) in Fig. El for Y = 8 and Y = 12, and in Fig. I14h for Y = 16. For Y = 8 where the solution is 
still symmetric, the difference between f(y, k 2 ) (lines) and / BK (y, k 2 ) (points) is quite large and low k 2 and 
visible at large k 2 , see Fig. 112b . The difference is significantly reduced at Y = 12, as clearly seen in Fig. 112b. 
Here, f(y, k 2 ) and / BK (y, k 2 ) almost exactly coincide except of some deviations for very small k < 0.1 GeV 
and y > Y/2. The overlap between f(y, k 2 ) and / BK (y, k 2 ) is further improved at Y = 16. Evolution of the 
smaller component /^ BK (y', k 2 ) in the BK limit may be also performed by solving the system (|27|) and (|28JI 
with the terms neglected that were generated by the triple pomeron vertex corresponding to the pomeron 
splitting (the contribution to the action of C\). The comparison of f^(y',k 2 ) and /^ BK (y',/c 2 ) is given for 
Y = 8, Y = 12 in Fig. and for Y = 16 in Fig. 114b . In this case, the two different kinds of solutions 
coincide even better than the large components f(y,k 2 ) and / BK (y,/c 2 ). 

Recall, that we also observed the similarity between the asymmetrical classical solutions of RFT-0 and 
the solution to the "fan equation", the counterpart of the BK equation in zero transverse dimensions. Thus, 
the "fan dominance" at Y S> Y c seems to be a generic feature of the interacting pomeron system. Even 
more can be said - the dependence of the QCD pomeron fields on the momentum might even enhance 
the convergence to the fan dominated system, compare Fig. and Fig. 1141 It happens probably because 
the deviations from the fan behavior appear at rapidities y — > Y, close to the source of the smaller field 
(which we chose to be P(y,k 2 ) for QCD pomerons and p(y) for RFT-0) where the smaller field is not yet 
strongly suppressed by the evolution. In QCD, however, the input is localized at rather small values of gluon 
momenta k, whereas the larger field, f(y,k 2 ), is concentrated around the saturation scale Q s (y) which is 
large for y —> Y. Therefore, the relatively large f^(y', k 2 ) in this rapidity domain affects only the tail of low 
momenta in f(y, k 2 ), with little relevance for the dynamics of the system. Possible implications of the "fan 
dominance" are discussed in the Sec. 03 

As the last point, let us comment shortly on the issue of multiple solutions to the Braun equations 
out of which only some can be found by the iterative solving procedure. Recall, that in the semi-classical 
approximation of the system moving in the Euclidean time, the most relevant are the trajectories with the 
lowest value of the action. We have no proof that the asymmetric classical trajectories that were found 
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in this paper fulfill this requirement. Numerous similarities of the patterns of solutions obtained in the 
interacting QCD Pomeron Field Theory and RFT-0 are, fortunately, reassuring. In both cases there exists 
a symmetric solution at low rapidity and two asymmetric solutions at Y > Y c . In both theories the "fan 
dominance" phenomenon was found. Therefore, one may conjecture that the asymmetric solutions of the 
Braun equations, indeed, are the classical trajectories with the lowest action, in analogy to the explicit result 
obtained in RFT-0. 

4.4 Summary of the results 

Let us summarize the presentation of results with a recapitulation of the most important observations: 

1. Solutions {/, /t} of the Braun equations with symmetric boundary conditions split into two different 
types: the symmetric solution f(y, k 2 ) = fi(Y — y, k 2 ) that dominates below the critical rapidity Y c 
and a pair of the asymmetric solutions, found for Y >Y C . 

2. The asymmetric solutions exhibit the feature of "fan dominance" which becomes more accurate with 
increasing Y. Due to the smallness of one of the field the system evolves as if one of the triple pomeron 
vertices (describing splitting or merging) was absent. The larger field is close to the solution of the 
BK equation. 

3. Unitarity corrections for the smaller field f\y\k 2 ) are very pronounced at Y 3> Y c leading to very 
strong damping (even 2-3 orders of magnitude at Y = 16, and increasing with Y) of the smaller field 
and flattening of the shape of fi(y', k 2 )/k 2 . We find that ft(y', k 2 ) ~ k 2 for k 2 < Q 2 (y), where Q s (y) 
is the saturation scale generated by the larger field f(y, k 2 ). 

4. The symmetric solution below the critical rapidity is significantly flatter at low momenta than the 
BK solution. At large momenta the decrease of the symmetric solution is power-like, with an exponent 
close to that of BK, but the solution to the Braun system is somewhat smaller the the BK solution 
from the same input. 

5 Discussion 

The breaking of the projectile-target symmetry which we have found above the critical rapidity is rather 
surprising and it calls for an explanation and interpretation. To our understanding the mechanism of this 
breaking is the following. Suppose that we have a symmetric situation in the system of the two evolving 
pomeron fields {/,/'}• If the fields are small and therefore weakly interacting then the interaction is only 
a small perturbation and the symmetry of the action and of the initial conditions should be reflected in 
the solution. This is, indeed, the case for the Braun equations with the total rapidity Y smaller than the 
critical value Y c . Let us consider now a symmetric system of pomeron fields when the fields are already 
strong due to their rapidity evolution. Then, the fields absorb intensively each other. The combination of 
the multiplication of the fields and the strong mutual absorption is a potential source of instability. Namely, 
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if we perturb the symmetric system of fields in this regime by, say, small increase of the value of the field / 
then the absorption of by interaction with / will be also increased. Thus, after this perturbation p 
should become smaller. This results, however, in a smaller absorption of the field /, leading to yet higher 
values of /, so that the instability will self-amplify generating finally an asymmetric configuration. Of course, 
whether this scenario is realized, depends on the particular form of the action. From our numerical results 
we conclude that this is, indeed, the case for the interacting pomeron fields above the critical rapidity in 
the classical approximation. It is curious that in the earlier study of the Braun equations the symmetry 
breaking was not found [HE]- Instead, there was reported an instability of the iterative procedure at critical 
values of rapidity, depending on the input. It was interpreted as a possible indication of a phase transition. 
We speculate that those instabilities might be, in fact, signs of emergence of asymmetrical solutions. 

The key question arises what happens with the spontaneous symmetry breaking (SSB) when the quantum 
effects are considered. Strictly speaking, for finite systems SSB does not occur at the quantum level. The 
ground state of the finite system in which the symmetry is spontaneously broken at the classical level is a 
symmetric coherent superposition of non-symmetric states. It is only in the thermodynamic limit when the 
SSB may take place in a quantum system. In the case of the reggeon field theory situation is even more 
complicated due to the fact that the evolution variable (the rapidity) would correspond to imaginary time 
in the Schrodinger picture. The consequences of the fact were investigated in detail in the framework of 
the reggeon field theory with impact parameter dependence (RFT-b) 47 . It turns out that in RFT-b the 
degenerate vacua communicate by quantum evolution irrespectively to the extension of the system in the 
transverse space and the symmetry of quantum theory is maintained even in the thermodynamical limit. The 
communication was found to be realized by solitons in the impact parameter plane smoothly interpolating 
between the two asymmetric vacua. In addition, the bifurcation of the classical solution at Y = Y c would 
indicate presence of a singularity of the S'-matrix 3 in Y. This singularity is not expected to be present in 
the complete quantum theory. 

It should be stressed, however, that in this paper we do not address the issue of properties of the QCD 
pomeron field theory in the thermodynamical context. The goal is rather to get insight into scattering of two 
strong sources of colour field e.g. the nuclei. In collisions of two nuclei the measurements give access not only 
to the total cross sections but also to extended information about the kinematics of the produced particles on 
the event-by-event basis. This is a classical measurement which, necessarily, breaks the quantum coherence. 
Therefore it is possible that such measurement selects one of the classical pomeron field trajectories which 
exhibit the symmetry breaking between the target and the projectile. In fact, the rapidity dependence of the 
saturation scale is different for the two asymmetric solutions of the Braun equations: for one of the solutions 
the saturation scale increases from the target to the projectile, while for the other it decreases. The average 
transverse momentum pT of the emitted particles should be correlated with the saturation scale. Hence, a 
classical measurement of the event should select one of the asymmetric solutions and it could exhibit some 
asymmetry in rapidity distribution of the produced particles. The pattern may be somewhat obscured, 
however, when the dependence on the transverse position is taken into account. In the collision the regions 
separated in the impact parameter are only weakly correlated and, in principle, it is possible that different 
3 We thank Lev Lipatov for this point. 



21 



domains in the transverse space are dominated by different asymmetric solutions of Braun equations. This 
would make the possible effects of asymmetry more subtle and harder to disentangle. 

One of the question which should be addressed is what observables could serve as experimental signatures 
of the asymmetry between the target and the projectile in heavy ion collisions. Certainly, the total cross 
section carries no information about the details of the evolution, so one should focus on more detailed 
observables. As a first guess we would propose investigation of the average transverse momentum pT = 



Y {p\) of the particles produced in central collisions of heavy ions as a function of rapidity y in the c.m.s. 

frame on the event-by-event basis. With the symmetry between the target and the projectile being preserved 
the observable Pt(u) measured for individual events should be the same after changing the definition of 
rapidity y — > —y. If the symmetry is broken in the event, however, Pt(u) should exhibit a clear trend. 

At this stage, we are not able to determine whether the symmetry breaking is a real physical phenomenon 
or an artifact of the effective theory of interacting pomerons. Needless to say, the framework of Braun 
equations relies on several assumptions that are far from being proven. First of all, it is not clear if the 
pomerons are valid degrees of freedom in dense and strongly interacting gluonic systems. One may argue 
that at high density the pomerons overlap and melt down to gluons, whose dynamics may be significantly 
different from the dynamics of the pomeron fields. Secondly, in the present analysis we neglected quantum 
effects related to the pomeron loops. The impact of the quantum effects on the phenomenon of symmetry 
breaking is unknown. Moreover, we neglected contribution of vertices with more than three pomerons. In 
addition, the NLL corrections to the BFKL pomeron kernels and to the triple pomeron vertices are neglected 
in the present form of Braun equations. This causes the BFKL intercept to be roughly two times too large. 
This means that the spontaneous breaking of the projectile-target symmetry should occur (if it occurs) at 
much higher rapidities than it may be deduced from the analysis employing the LL BFKL kernel, perhaps 
for energies beyond reach of the LHC. On the other hand, the value of a s = 0.2 underestimates significantly 
the expected value of a s for the triple pomeron vertex (recall that the vertex is proportional to a 2 .) and with 
a more realistic value a smaller Y c would be predicted. Finally, in order to evaluate relevance of the effect 
the input conditions should be carefully tuned to embody the available information on the unintegrated 
gluon density in the nuclei, including the impact parameter profiles. Keeping in mind all these reservation, 
we believe that further theoretical and experimental studies of the issue should be carried out. 

Leaving the issue of the symmetry breaking, we point out that the "fan dominance" phenomenon at 

Y > Y c found in the case of the symmetric boundary conditions should be even more pronounced when 
the projectile and the target are different. This might provide some basis for the use of the BK equation 
to describe the saturation effects in the DIS at low Q 2 . Strictly speaking, the BK equation is valid for a 
small perturbative probe scattering off a large target, for instance a nucleus. This condition is, certainly, 
not fulfilled for almost real photon scattering off a proton, nor it is for the diffractive DIS, dominated by 
scattering of large dipoles. Still, the fits based on the BK amplitudes are very successful in both cases. The 
"fan dominance" in the symmetric Braun system could provide some support for those applications of the 
BK equation, although it is fair to admit that the use of Braun equations to processes of this kind is not on 
the firm ground either. 
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The Braun equations are a minimal extension of the very fruitful concept of the BK equation, that 
embodies the symmetry between the pomeron fields / and /' at the level of the action. Possibly, this 
extension may also find some interesting phenomenological applications. Following Braun we state that 
description of heavy ion collisions is the obvious application of the equations. In that case the quantum 
loops should have only a subleading effect and the approximate treatment of the dependence transverse 
position is certainly sufficient. A more challenging is an application of the formalism to the vital problem 
of understanding of pp collisions at the LHC energies. In particular, we have in mind the description of 
underlying event, particle production (see e.g. [Ml EH])) diffractive processes and determination of the gap 
survival factor in hard exclusive processes, like for instance the exclusive Higgs boson production. Expanding 
on this example, the exclusive Higgs boson production is an important process which is probable to be 
measured at the LHC It was shown, that the theoretical understanding of the cross section for this 

process requires a good control of the hard rescattering corrections |51| . The framework of the semi-classical 
field theory of the interacting pomerons may serve as a tool to perform the necessary resummations of multi- 
pomeron diagrams and to obtain improved estimates of the gap survival factor. All the listed applications 
are, however, non-trivial as they require inclusion of NLL BFKL corrections and, possibly, more accurate 
treatment of the dynamics of the system in the transverse position space. 

In the last part of this section we will briefly mention some intriguing open questions. Thus, it would be 
interesting to investigate the stochastic QCD evolution of the color glass condensate using the realization of 
the semi-classical approximations in which only the tree topologies of the pomeron are retained. In analogy 
to the Pomeron Field Theory, this limit should be simpler than the accurate treatment also in the CGC 
formulation. Furthermore, a similar analysis of the Pomeron Field Theory should be possible after inclusion 
of the pomeron vertices at which more than three pomeron fields. The form of those vertices may be 
predicted using the conjectured conformal symmetry of the pomeron field theory j52j . While the conformal 
symmetry of the effective field theory of the interacting pomerons in the EGLLA was not explored in the 
present study, it constitutes, certainly, a key ingredient of the structure of the complete theory. Thus, it is 
mandatory to account for it in similar future studies. 

Finally, let us refer to recent developments on the connection between the superconformal gauge theories 
in four dimensions and the superstring theory on the AdS^ x S 5 background pQ. Using the AdS/CFT 
duality it was found that the BFKL pomeron in gauge theories corresponds to the graviton Regge trajectory 
in the AdS space [HH]- Thus, it is desirable to find an interpretation of the effective pomeron field theory at 
the string side, perhaps in terms of gravity. Curiously enough, it was discovered recently that collisions of 
heavy ions possess a dual gravitational description jSl]. The dual of the scattering is given by a collision of 
two gravitational shock waves in which black holes can be formed. Example of such a black hole solution 
being produced, that moves in the fifth dimension of the Anti de Sitter space was found [53] • Thus, the 
relation between the fifth dimension in the AdS and the gluon virtuality inspires a question about the 
possible connection of the black hole to the BK traveling wave solution, whose existence we established in 
the classical pomeron field theory. Therefore, it is important to verify whether phenomena analogous to 
the symmetry breaking between the target and the projectile and the "fan dominance" also happen in the 
string world. 
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6 Conclusions 



Effective field theory of interacting QCD pomerons was investigated in the semi-classical limit, as a frame- 
work to describe high energy scattering of two nuclei. The effective action was proposed in the form using 
the pomeron amplitudes, as the basic degrees of freedom, related to the unintegrated gluon densities in the 
linear regime. Triple pomeron vertices in the momentum space accounted for the pomeron merging and 
splitting. Arbitrariness in the choice of the direction of evolution in rapidity required both the vertices to 
be identical and induced the self-duality of the action. This symmetry combined with symmetric initial 
conditions defined the scattering problem to be symmetric under the interchange of the target and the 
projectile. Classical pomeron field equations {Brawn equations) were re-derived and solved numerically. The 
solutions were found that were invariant under the projectile-target symmetry only for scattering rapidi- 
ties Y smaller than a critical (non- universal) value Y c . For Y > Y c the projectile-target symmetry turned 
out to be spontaneously broken. Above the critical rapidity, the solutions converged to the solutions of the 
Balitsky-Kovchegov equation, the phenomenon which we called the BK fan dominance. A very similar pat- 
tern of symmetry breaking and the fan dominance occurs also in the Reggeon Field Theory in zero transverse 
dimensions, which suggests that this is a generic feature of the interacting pomeron system. We discussed 
possible consequences of those observations for the phenomenology of heavy ion collisions and the physics 
of pp scattering at the LHC. Finally, we suggested that the results of this paper may have counterparts in 
the dual description of heavy ion collisions in terms of scattering of two gravitational shock waves in the 
Anti de Sitter space in five dimensions. 
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Figure 7: Solutions of the Braun equations f(y,k 2 )/k 2 for a) Y = 10 and b) Y 
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Figure 8: Solutions of the Braun equations f* (y', k 2 )/k 2 for Y = 10: a) y' = 0, 1, . . . , 5; b) y' = 6,7, . . . , 10. 
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Figure 9: Solutions of the Braun equations f^(y', k 2 )/k 2 for Y = 16: a) y' = 0, 1, . . . , 8; b) y' = 9, 10, . . . , 16. 
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Figure 10: Solutions of the Braun equations for Y = 16 plotted as a function of rapidity y for k = 1.3 GeV 
(solid line), k = 14 GeV (dashed line) and k = 185 GeV (dotted line): a) f(y,k 2 )/k 2 and b) P{y, k 2 )/k 2 . 
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Figure 11: Solutions of the Braun equations for Y = 16 plotted as a function of rapidity y and k: 
a)f(y,k 2 )/k 2 and b) p(y,k 2 )/k 2 . 
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Figure 12: Comparison of the larger solution f(y,k 2 )/k 2 of the Braun equations (solid line) to the solution 
of the BK equation (points) for a) Y = 8 and b) Y = 12. 
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Figure 13: Comparison of the smaller solution f^(y',k 2 )/k 2 of the Braun equations (solid line) to the 
"smaller solution" f' (BK (y',k 2 )/k 2 of the BK equation (points) for a) Y = 8 and b) Y = 12. 
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Figure 14: Comparison of the solution of the Braun equations (solid line) to the solution of the BK equation 
(points) forY = 16: a) the larger solution f(y,k 2 )/k 2 and b) the smaller solution p(y',k 2 )/k 2 . 
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